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Abstract 

A model potential for two-particle relativistic systems is investigated in the 
framework of Poincare-invariant quantum mechanics (or relativistic Hamiltonian 
dynamics). The potential considered allows to reduce the main 
integro-differential equation of Poincare-invariant quantum mechanics to the 

equation analogous to the radial equation and have analytical solution for 
relativistic bound system. We discuss the possible choice of the parameters of 
potential and apply our model to the description of light meson Regge 

trajectories. 
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1 Introduction 

Relativistic few body problem has received a great 
attention in hadronic and nuclear physics. The 
most complete results here exist for the case of two 
particles. Description of the bound system in the 
relativistic quantum field theory is founded on the 
four-dimensional covariant Bethe-Salpeter equation 
However, this equation gives series of difficulties 
when the practical calculations are made. 

There exist various reductions of the two-body 
Bethe-Salpeter equation. Different forms of this re- 
duction were discussed Logunov and Tavkhelidze 
§, Kadyshevsky §, Todorov @, Gross §, Polu- 
zou, Keister and Lev and many others. 

The relativistic two-body systems were analyzed 
by Weinberg Q, Frankfurt and Strikman |K|, Kon- 
dratyuk and Terent'ev (ll]] in infinite-momentum 
frame. Some authors use diagrammatic approach 
i.e. they select leading diagrams and project them 
onto the three-dimensional space. Others make use 
effective Hamiltonians, Faddev equation. Different 
forms of the quasi-potential equation can be derived 
using this approach 

In this work we apply the point form of rela- 
tivistic quantum mechanics (RQM) @-|8[ to the 
description of meson Regge trajectories. RQM is 



also known in the literature as relativistic Hamil- 
tonian dynamics or Poincare-invariant a quantum 
mechanics with direct interaction. 



2 RQM formalism for quark- 
ant iquark bound states 

The formulation of relativistic quantum mechanics 
differs from non-relativistic quantum mechanics by 
the replacement of invariance under Galilean trans- 
formations with invariance under Poincare transfor- 
mations. The dynamics of many-particle system in 
the RQM is specified by expressing the ten genera- 
tors of the Poincare group Muu an d I « in terms of 
dynamical variables. In the constructing generators 
for interacting systems it is customary to start with 
the generators of the corresponding non-interacting 
system (we shall write this operators without "a 
hat" ) and then add interaction in the a way that is 
consistent with Poincare algebra. In the relativistic 
case it is necessary to add an interaction V to more 
than one generator in order to satisfy the commu- 
tation relations of the Poincare algebra. Dirac [14| 
observed, that there is no unique way of separat- 



ing the generators into dynamical subset (the gen- 
erators including interaction V) and kinematic sub- 
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set. Kinematic subset must be associated with some 
subgroup of Poincare group, usually called stability 
group or kinematic subgroup. Thus, the con- 
struction for two interacting particles proceeds as 
follows: 

1. The two-particle Hilbert space of non- 
interacting system is defined as tensor product of 
two one-particle Hilbert spaces. A two-body uni- 
tary representation of Poincare group on the two- 
particle Hilbert space is reducible. A basis in this 
space can be constructed from single-particle bases: 



\pi\i) \p 2 X 2 ) = |misi;piAi) <8> \m 2 s 2 ;p 2 X 2 ) (1) 
with normalization 

(p'jXI <P 2 A 2 | IfiAi) |P2A 2 ) 

where Pi, mi, Sj, Aj are momenta, masses, spins and 
projects of spin of particles, which will form bound 
system. 

2. Clebsch-Gordan coefficients for Poincare group 
are constructed and used to reduce the unitary rep- 
resentation on the two-particle Hilbert space to lin- 
ear superposition (direct integral) of irreducible rep- 
resentations. Poincare generators for irreducible 
representations of the non-interacting system are 
constructed, along with operators for the mass and 
spin and other operators. The result of this step is 

P 12 , (i, [J, k] , (Is); [misi; m 2 s 2 ] 



EE/ d\ 



U mi (Pl)um 2 (P2)M 



is A1A2 " V u mi (k)u m2 (k)u Mo (Pi2) 
E E <si^i> s 2 ^2| sX) (lm,s\\ JfJ,)Y lm (9,. 



rn\ vwi 
.1/2 



Dj£ (n( Vl ,P l2 ))D l H 2 (n(p 2 ,P 12 )) \p1X1) \p 2 X 2 ) , 

(2) 

where (s\V\, s 2 i/ 2 | sX), {lm, sX\ Jfj) are Clebsh- 
Gordan coefficients of £J7(2)-group, YJ m (#,</>) are 
the spherical harmonics with spherical angle of k. 
Also, in Eq. (|) D 1 / 2 (n) = 1 - iti-a/y/l + n* is D- 
function of Wigner rotation, which is determined by 



the vector-parameter n(pi,p 2 ) = u\ xu 2 /(l— u\ -u 2 ) 
with u = pi (uj m (p) + m). The three- momenta p±, 
p 2 of the particles (quarks in our case) with the 
masses mi and m 2 of relativistic system are trans- 
formed to the total, P12, and relative momenta k to 
facilitate the separation of the center-of-mass mo- 
tion: 

P12 = Pi +P2, 



k = pi + 



where 



P 



1.2 



P 



12 - Pi 



M \UM (P12)+M 



+ %1 (pi 



(3) 



M = u) mx (k) +0J m2 (k) (4) 

mass of non- interacting system and oJ mi (pi) = 

p[ + m\ . 
The 

state vector P12, /x, [J k] , (I s);[misi;m 2 s 2 ] 

P12, J, fJ-, k, (I s)^ is the eigenstate of operators P12 

, J 2 , J3 and, also, L 2 and S 2 , where L and S are rel- 
ative orbital momentum and total spin momentum, 
respectively. This vector of the non-interacting 
Qq system transforms irreducibly under Poincare 
transformations. 

3. Following Bakamjian and Tomas, interactions 
are added to the mass operator in the irreducible 
free-particle representation. The new mass operator 
is defined by 

M = M + V. (5) 

If V is any operator that satisfies the following con- 
ditions 

M = M\ M > 0, (6) 

[P12, V}- = \iVp 12 , V}- = [J, V}- = (7) 

then the similar set of interacting particles will sat- 
isfy the same commutation relations as the set of 
non-interacting system. 

There are three forms of the dynamics in the 
relativistic quantum mechanics called "instant", 
"point", and "light-front" forms Ijjlf] . The descrip- 
tion in the instant form implies that the opera- 
tors of three-momentum and angular momentum 

do not depend on interactions, i.e. P = P and J 
= J (J = (M 23 , M 31 , M 12 )) and interactions can 
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be presented in terms of operator P° and genera- 
tors of the Lorentz boosts N = (M 01 , M 02 , M 03 ). 
The description in the point form implies that the 
operators M^ v are the same as for non-interacting 
particles, i.e. M^ v = M^ v , and these interaction 
terms can be presented only in the form of the four- 
momentum operators P. In the front form with the 
fixed z axis we introduce the + and - components 
of the four- vectors as p + = (p° + p z )/\^2, p~ = 
(p° — p z )j\pl. We require that in the front form 
the operators P + , Pi , M 12 , M+~ , M+i (j = 1,2) 
are the same as the corresponding free operators, 
and interaction terms can be presented via the op- 
erators M~ 3 and P~ . 

In our work we use the point form of RQM be- 
cause: 

a). The meson wave function (see below) being 
Lorentz invariant, as (in a point form of dynamics) 
the operator of a boost does not contain interaction. 
In the instant form of RQM this property of wave 
functions does not take a place. 
(3). The relativistic impulse approximation in an 
instant form and dynamics on light front automat- 
ically breaks down the Poincare-invariance of mod- 
els. In the point form of RQM such a violation does 
not happen. 

Given such an interaction V it is useful to define 
two other related interactions Wi, 



U = M 



M 2 



V 2 + MqV + VM 



(8) 



W 



[M 2 - Ml) + (m( - ml) 2 * 



(. 



1 



Vm 2 m 2 ) 

(9) 

The eigenvalue problem for the mass of Qq system 
can be expressed in the three equivalent forms ffl]: 



+ 2\J 77(17 + m 2 + m|) + mfm 2 ,. (11) 
When mi = rri2 = m, the equation ([ll]) reduces 



to 



M 2 = 4 (77 + m' 



(12) 



Solution of any of the above eigenvalue problems 
in the point form of RQM leads to the eigenfunc- 
tions of the form 



(13) 



V 12 ,J^,k,(ls) V,J',n',M 
Sjj'5^5(V -V 12 )^(kls) 



with the velocities of the bound system V = P/M 
and non-interacting system V12 = Pu/Mq. The 
function \E , " r/i (k I s) satisfies (in the point form) the 
following equation (7|, 

00 

I < k (I s) || W J || k' (I' s') > ^ J {k' I 1 s')k' 2 dk'+ 

vs'{ 

k 2 ^ J (kls) =r]^ J (kls) (14) 
with reduced matrix element of operator W 



V l2 ,J,fi,k,(ls) W V{ 2 ,J'^',k',{l's') 



Sjj'S^5(V 12 -V{ 2 )(k, (ls)\ 



\k>,(l>s')). (15) 



Equation ( |l~4| ) is a radial equation. Therefore we 
have used Poincare group properties to separate the 
radial and angular dependencies. 

As vector of irreducible basis for a system with in- 
teraction and without its satisfies to conditions of a 
normalization and completeness, the wave functions 
satisfy to the following normalization condition, 



E 

Is 



dk k 



* J (kls) 



1. 



(16) 



M I f >= (M + V) I * >= M I * >, 

Ml + U I ^ >= M 2 I * >, 
{k 2 + W) I ^ >= rj I f >, 



(10) 



where the mass M of bound Qq system and r\ have 
relationship, 



M z 



2rj + m\ + m\+ 



3 Model potential and solving 
main equation. 

In this section we apply the formalism developed 
above to calculate mass spectra of mesons contain- 
ing u, d and s quarks. To choose appropriate in- 
terquark potential we use a well-known experimen- 
tal fact that light hadrons populate approximately 
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linear Regge trajectories, i.e., M 2 ~ (3 J + const, 
with the same slope, (3 — 1.2 GeV 2 , for all trajec- 
tories (see, for example, [[H]]). Therefore, we can 
take the effective model potential W in the oscilla- 
tor form 



has the Gaussian form, which is used in many rel- 
ativistic models of hadrons. This wave functions 
gives an excellent description of electromagnetic 
formfactors at small transfers Q 2 , Q 2 <a few GeV . 
Quantization condition is defined by 



W(r) = W 5(r) + f3 4 r 2 , 



(17) 



n = Wo + 2(3 2 (2n + l + 3/2) 



(20) 



where Wq and (5 are free parameters. As follows 
from calculations of meson spectra done by God- 



frey and Isgur [17], spin-dependent corrections are 
important for IS", 2S, l 3 i-2 and 1 3 P 2 states. How- 
ever for other states spin-dependent corrections are 
small and they can be neglected in Eq. (|l7|). Us- 
ing standard relationships for operator f and orbital 

momentum operator L 2 



(k\r\k') = -iV k 5(k-k'), 



V- 



V 12 ,J,fi,k, (Is)) = 1(1 + I) V 12 ,J,p,k, (Is) 



we reduce main integro-differential equation of 
RQM to the ordinary quantum-mechanical radial 
equation with the oscillator potential (only impulse 
representation), 



d 2 2d_ 

dk 2 k dk 



_ 1(1 + 1) 

k 2 

W -v 



/? 4 



V(kls) 



V 4 



ty(kls). 



(18) 



The eigenfunctions of Eq. (18) are 
= JV ni exp 




with 



nl 



-n.l + 



3 e 

2'/3 2 



(19) 



-1-3/2 



T(n + l + 3/2) 



r(z + 3/2) V r(n+l) 



where n, I = 0,1,2..., F(a,b,z) is the hypergeo- 
metric function, T(n) is the Gamma function. Note 
that the wave function of the ground state ( n, I = 0) 



The spectra of mesons composed of quarks with 
equal masses (mi = m-2 = m ), for example, u and 
d quarks, are given by (see Eq. (O) and Eq. (|20|) ) 



M 



4(m z + W ) + 8/3 2 2n + I + - . 



(21) 



Thus we reproduce the linear dependence M 2 (l) in 
the framework of the two-body relativistic equation 



(14). The parameters (m 2 + Wq) and (3 have been 
found from the fitting the Regge trajectories (see 
Fig.l), 
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Figure 1: p — meson Regge trajectory 

(m 2 + W ) = -0.261 GeV 2 , (3 = 0.375 GeV 
If we take the same masses of u and d quark as in 



Ref. i.e., 

m u = md = m = 0.22 GeV, 
then we obtain 

Wq = -0.31 GeV 2 . 



(22) 
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There are eight meson Regge trajectories popu- 
lated by u — d bound states (for each isospin / and 
angular momenta J = I rb 1, J = I and total spin 5 
of qq system, 5 = 0, 1. Some of the trajectories are 
plotted in Fig. 2-4. 

We observe that all experimental data are in good 
agreement with spectrum given by Eq. (|2l"|) for 
I > 1 and 5 = 1 (Fig. 2) . As in the case of bound 
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Figure 2: w-meson Regge trajectory 

states with 5 = 0, the agreement between our the- 
oretical predictions and the existing experimental 
data is not good (see Fig. 3-4). Such deviations can 
be explained by absence of spin-dependent terms 
and short-distance term of the potential. 

As it is easy to see, our method for bound systems 
with different masses of quarks does not require spe- 
cial processing of solving main equation of RQM, as 
it was made in Ref.|l(| (introduction of additional 
parameter). When the masses of the quark and an- 
tiquark are different we obtain 

M 2 = 2W + 2/3 2 (2n + I + 3/2) + m\ + m 2 2 + 



+2^{W Q + 2(3 2 (2n + l + 3/2))* 

* \/(W + 2j3 2 (2n + l + 3/2) +m\ + m\) + m\m\ 

(23) 

We see that the dependence (|23|) is, also, linear, 
but asymptotically at large I. We apply (23) to the 
description of strange meson trajectories. For the 
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Figure 3: 77-meson Regge trajectory 

mass of the u(d)— quark and for mass strange quark 
we used the same value as in Ref. Q, i.e., mi = 
m u = 0.22 GeV (see ) , rn 2 = m s = 0.419 GeV. 
As in the case of meson with hidden strangeness, 
the agreement between our model predictions and 



existing experimental data [18| is good (see Fig. 5). 



Also we can predict all the strange meson tra- 
jectories (see Fig. 6-9). The model describes all 
s — q(q — s) meson Regge trajectories with 5=1 
in quite satisfactory way (see Fig. 6-8). But for 
mesons with 5 = 0, the agreement between model 
and the experimental data is not good (see Fig. 9), 
as in the case of u — d mesons. 

In summary we shall mark that, the analyti- 
cal solution of a main equation a RQM (|l9| ) can 
be used as zero approximation for problem solv- 
ing with more realistic potentials. The model ad- 
equately describes, with only four free parameters, 
m u = rrid, rn s , Wo, (3, orbitally excited meson Regge 
trajectories. 
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